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We discuss entanglement evolution of a multi-qubit system when one of its qubits is subjected to 
a general noisy channel. For such a system, an evolution equation of entanglement for a lower bound 
for multi-qubit concurrence is derived. Using this evolution equation, the entanglement dynamics 
of an initially mixed three-qubit state composed of a GHZ and a W state is analyzed if one of the 
qubits is affected by a phase, an amplitude or a generalized amplitude damping channel. 
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I. INTRODUCTION 

Entanglement is an exclusive nonclassical resource of 
quantum mechanics with promising applications in com- 
munication, cryptography and computing In many 
of these applications, this resource is often required to 
be distributed between different partners; a request that 
implies that at least one physical (sub-)system is to be 
transmitted through a communication channel. In gen- 
eral, such a coupling of a quantum system to some en- 
vironmental channel leads to decoherence of the system 
and usually to some loss of entanglement before the state 
of the system can be further utilized. Therefore, it is of 
great practical importance to investigate the decoherence 
and time evolution of entanglement for quantum systems 
that undergo the action of noisy channels. 

Typically, the time evolution of entanglement of a sys- 
tem is deduced from studyingits state evolution under 
the influence of decoherence 043- Following this line, 
we have recently analyzed for example the entanglement 
dynamics of initially prepared pure three-qubit GHZ and 
W states if transmitted through one of the Pauli channels 
02! cTxj &y or a depolarizing channel [BJ . In this work, we 
employed the analytical solutions of the master equation 
as obtained by Jung et al. IfJ and a lower bound for 
the three-qubit concurrence [7| in order to quantify the 
time-dependent entanglement of the mixed three-qubit 
states. 

Instead of making explicit use of the state evolution 
for the analysis of the entanglement dynamics of a given 
system, Konrad et al. [|| recently derived an evolution 
equation of entanglement for a two-qubit system that pro- 
vides a direct relationship between the initial and the fi- 
nal entanglement of the system when one of its qubits 
is subjected to an arbitrary noise. Subsequently, Li et 
al. |9| derived a generalized evolution equation of en- 
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tanglement for a (finite dimensional) bipartite system, if 
initially prepared in a pure state and affected by an arbi- 
trary noisy channel. Recently, moreover, this latter result 
Q has been extended to the case of an initial mixed state 
of a bipartite system fiol . [Tl| . 

In this work, we suggest an evolution equation of en- 
tanglement for a multi-qubit system when one of its 
qubits undergoes the action of an arbitrary channel, 
which is given by a completely positive (non-)trace- 
preserving map. Since there is no an analytically com- 
putable measure of entanglement for mixed multi-qubit 
states [lH , we shall make use of the lower bound by Li et 
al. Q for the multi-qubit concurrence. In fact, this lower 
bound was proven to satisfy the general requirements Q 
of an entanglement measure 0, [Hf ■ Therefore, the evo- 
lution equation derived below is capable to characterize 
the entanglement dynamics of multi-qubit states repre- 
senting a lower bound for an actual currently unavailable 
entanglement evolution. As example, we shall discuss in 
details the time evolution of entanglement of an initially 
mixed three-qubit state composed of a GHZ and a W 
state Q 

p(p) = p \GHZ) (GHZ\ + (1 - p) \W) (W\ , (1) 

if one of the qubits is affected by a phase or an amplitude 
damping channel. For a generalized amplitude damping 
channel, moreover, we also show the sudden death of en- 
tanglement [I| for this mixed three-qubit state. 

The paper is organized as follows. In the next section, 
we first recall how the entanglement of a multi-qubit sys- 
tem can be quantified either by means of the iV-qubit 
concurrence 0, HH or in terms of a lower bound to this 
measure as suggested by Li et al. Q. This lower bound 
was constructed in such a way that it includes only bipar- 
tite concurrences according to some "bi-partite" cuts of 
the multi-qubit system. In Section. Ill Bl we display and 
discuss the evolution equation of bipartite concurrence as 
obtained by Li et al. @ . Based on their results for the bi- 
partite concurrence, we then derive an evolution equation 
of the lower bound to the concurrence for a three-qubit 



2 



system in Section II. C, if one of its qubits is affected by an 
arbitrary noisy channel and if we start with an initially 
pure state. In Section Hi Dl we shall discuss also possi- 
ble extensions of this equation to the cases of N qubits, 
initially mixed states of three qubits, and if a system is 
affected by many-sided noisy channels, i.e. when several 
of its qubits undergo simultaneously the action of some 
local noise. In Section IIII1 we later analyze in detail the 
entanglement dynamics of the initially thrcc-qubit mixed 
state (JT|) for three realistic noise models, namely, a phase 
damping, an amplitude damping and a generalized am- 
plitude damping noise. Finally, a conclusion is drawn in 
Section [TV] 



II. EVOLUTION EQUATION OF THE LOWER 
BOUND FOR MULTI-QUBIT CONCURRENCE 

Until the present, it has been found difficult to quan- 
tify the entanglement of mixed many-partite states, and 
no general solution is known |l2| apart from Wootter's 
(two-qubit) concurrence [HI]. This concurrence provides 
indeed a very powerful measure of entanglement but is 
just suitable for two-qubit systems. Various extensions 
of Wootter's concurrence have been worked out over the 
years, and especially for mixed bipartite states, if the di- 
mensions of the associated Hilbcrt (sub-)spaces are larger 
than two [H, [III G3 • Nonetheless, a full generalization of 
the concurrence towards mixed many-partite states has 
remained a challenge until now and will require further 
studies [l^ . In practice, however, one is mostly inter- 
ested in the minimum (nontrivial) amount of entangle- 
ment fl8j which is preserved in a mixed state of some 
system when coupled to its environment, for instance, 
due to transmission of one or several of its subsystems 
through a communication channel. To this end, vari- 
ous analytically computable lower bounds to concurrence 
were suggested recently 0, 0, EH . In the next section, we 
shall exploit the lower bound to the multi-qubit concur- 
rence as suggested by Li et al. Q. As mentioned before, 
this bound is based on the bipartite concurrences of the 
multi-qubit system with regard to some bipartite cuts. 
This "bi-partitioning" enables us eventually to construct 
an evolution equation for multi-qubit systems which is 
based on the evolution equation for bipartite systems Q . 



A. A lower bound for iV-qubit concurrence 

For a given pure N-qubit state the concurrence 
can be written as @, H(| 



the remaining N — 1 qubits. Since any mixed state can 
be expressed also as a convex sum of some pure states 
P = J2i Pi IV'i) the definition (J2j for the con- 
currence for pure N-qubit states can be generalized for 
mixed states by means of the so-called convex roof [l2j 



C N {p) = min p.j C 



(IV>i» 



(3) 



In this latter expression, however, the minimum has to 
be found with regard to all possible decompositions of p 
into pure states \ipi). 

Unfortunately, no solution has been found so far to 
optimize the concurrence ([3]) of a multi-qubit system an- 
alytically [HI , apart from the two-qubit systems [lj| and 
a special case of three qubits [l4[ . Instead of the compu- 
tationally demanding optimization of the right-hand side 
of Eq. ([3]), a much simpler lower bound to this measure 
has been suggested recently by Li et al. Q 



Cjv(p) > Tjv(p) 
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N K 



(4) 



This bound is defined in terms of the N "bi-partite" con- 
currences C n [HI that correspond to the possible (bipar- 
tite) cuts of the multi-qubit system in which just one of 
the qubits is discriminated from the other N — 1 qubits. 
For the separation of the n-th qubit, the bipartite con- 
currence C n is given by a sum of K — 2 N ~ 2 (2 iv_1 — 1) 
terms C& which are expressed as 



C£=max{0,A£-A 2 



(5) 



and where the A 



fe i 



1..4 are the square roots of 
the four nonvanishing eigenvalues of the matrix p p 1 ^ , if 
taken in decreasing order. These (non-hcrmitian) ma- 
trices p pj? are formed by means of the density matrix 
p and its complex conjugate p* , and are further trans- 



formed by the operators {S% 



1, 



,K} 



as: p% = Sup*S^. In this notation, moreover, Lo is the 
(single) generator of the group SO(2), while the {-£<£} 
are the K = 2 N ~ 2 (2 iv_1 — 1) generators of the group 
SO(2 Af_1 ). We note that the lower bound (0} reduces 
to the Wootters concurrence [l5| for just two qubits. For 
details about the explicit construction of the lower bound 
dH) we refer to 0- 

Let us display this lower bound (0| especially for three- 
qubits, Tz(p), for which we wish later to describe the en- 
tanglement dynamics. For such states, the lower bound 
T${p) can be written in terms of the three bipartite con- 
currences that correspond to possible cuts of the two 
qubits from the remaining one, i.e. 
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(cf 13 ) 2 



(cf 2 ) 2 



(cf 1 ) 2 , 



(G) 



and where the pi = Tr \ip) (ip\ denote the reduced density 
matrix of the i-th qubit which is obtained by tracing out 



The bipartite concurrence C^ b ' c (for a, 6, c = 1..3 and 
a =/= b =/= c =/= a) arc obtained as described above with the 
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help of the operators {S^ c = L^^Lq, k = 1...6}, where 
L is the generator of the group SO (2) which is given by 
the second Pauli matrix a y = -i (|0) (1| + |1) (0|). The 
(six) generators L^ b of the group SO (4) that can be ex- 
pressed explicitly by means of the totally antisymmet- 
ric Levi-Cevita symbol in four dimensions as (Lki) mn = 
—i £kimn] k, l,m,n = 1..4 [HJ. 

B. Evolution equation for bipartite systems 



|0) of the bipartite system if the singlc-qubit subsystem 
is affected by the noisy channel S. 

We can utilize Eq. ([9]) to derive next an evolution equa- 
tion of the lower bound ([5]) to the threc-qubit concur- 
rence, if just one qubit is affected by some noisy channel. 
Later, wc shall generalize this evolution equation to the 
case of a general ./V-qubit system with the same assump- 
tion that just one of its qubits is subjected to a noisy 
channel. 



Having the explicit expression (0| for the lower bound 
tn(p) to the multi-qubit concurrence, we can derive an 
evolution equation of this entanglement measure. Since 
the lower bound Tjy(p) is given just in terms of the bipar- 
tite concurrences C n , let us first reconsider the evolution 
equation for a bipartite system as suggested by Li et al. 
[9[. Suppose, |x) is a pure state of a bipartite system 
d\ ® di with dimension d\ and d2 of the corresponding 
subsystems, and the second subsystem undergoes the ac- 
tion of a noisy channel S. Then, the final state of the 
system is a mixed state in general and takes the form 
p = (1 ® <S) |x) (x|. On the other hand, any pure state 
\x) can be obtained also from the maximally entangled 

state |0) = Yli=i Ijl® \i) /Vd~2 of the bipartite system by 
|x) = (M® 1) |0) Q|. I n this notation, M denotes a local 
(filtering) operator |22l . |23| that acts on the first subsys- 
tem of the maximally entangled state. Although local 
operations, usually associated with projective measure- 
ments or a passage through a channel, cannot change en- 
tanglement because of its nonlocal nature, there is a spe- 
cial case of a stochastic local operation, a filtering opera- 
tion, that can be employed to influence on this nonlocal 
feature [23| . Therefore, the final state p of the bipartite 
system can be expressed as 

p = (1 ® S) (M ® 1) |0) (0| (Mf ® 1) . (7) 

Since the filtering operator M and the noise S act only 
on either the first or second subsystem, it has been shown 
in Q that the reduction of the entanglement of the sys- 
tem under the action of a noisy channel S is independent 
of the initial state |x), and is bounded from above by the 
channel's action upon the maximal entangled state |0). 
We therefore obtain 

C[(l ® S) |x) (X|] < y C[(l ® S) |0) (0|]C[|x)] , (8) 

where C [..] denotes the bipartite concurrence 
If, moreover, the bipartite system consists of a tri- 
dimensional and a single-qubit subsystem, and just the 
qubit is affected by the noisy channel iS, the equal sign 
applies in inequality (|5J and wc obtain 

C[(l ® S) |x) (x|j = C[(l ® S) |0) (0|]C[|x)] . (9) 

That is, the entanglement dynamics of an arbitrary pure 
state of a d\<S>2 bipartite system is completely determined 
by the channel's action on the maximally entangled state 



C. Dynamics of the lower bound for initially pure 
three-qubit states 

Suppose |x) is a pure state of a three-qubit system 
and just one qubit undergoes the action of a channel S. 
The final state of the three-qubit system takes the form 
p = (1®1®<S) |x) (xIj which is equivalent to the final state 
of a bipartite 4 ® 2 system when the second subsystem 
is subjected to the channel S. As we mentioned above, 
any pure state of a bipartite system can be obtained from 
the maximally entangled state of the bipartite system by 
means of a single local filtering operation M acting on 
the first subsystem as |x) = (M® 1) |0). In contrast, two 
local filters M and M' are in general required to obtain 
an arbitrary pure three-qubit state |x) from a maximally 
entangled state of three qubits |0) by |x) = (M ® M' ® 
1) |0). For three-qubit systems, however, there are two 
maximally entangled states which, in the computational 
basis, are given by [24| 

|GHZ) = -^(|000> + |111» , (10) 
\W) = -J= (|001) + |010) + |100)) . (11) 

These two entangled states cannot be obtained from each 
other by means of local singlc-qubit unitary operations 
[24] and give rise to two local unitary inequivalent classes 
of three-qubit entangled states, the so-called GHZ- and 
W-classes. 

Although an arbitrary pure three-qubit state |x) can 
be generated from one of the maximally entangled states 
by means of local operations [24| , we first need to iden- 
tify the class of states either (|10p and pT|) to which it 
belongs to. For an arbitrary (pure or mixed) three-qubit 
entangled state, fortunately, this is possible by follow- 
ing the procedure due to Diir et al. [24j which is sim- 
ple and just includes the computation of the 3-tanglc as 
described in Ref. [l6j]. It leads to the distinction that 
every entangled three-qubit state |x), for which the 3- 
tangle vanishes, belong to the W-class and can thus be 
obtained from the W state (fTTj) by means of local uni- 
tary operations. In contrast, any entangled three-qubit 
state with nonvanishing 3-tangle is part of the GHZ-class. 
For a given pure three-qubit state |x), it is therefore al- 
ways possible to find proper local (filtering) operations 
M and M' so that |x) is obtained from either (|TU1) or 
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(fTTj) by |x) = (M M' 1) Moreover, we have 
|0) = | GHZ) if |x) belongs to the GHZ-class of entangle- 
ment, and \(f>) = \W) for |x) being part of the W-class. 

To summarize our discussion here, the final state of the 
three-qubit system when one of its qubits undergoes the 
action of a noisy channel S is given by 

p = (1 ® 1 ®<S) x 

{M®M' ® 1) |0) (0| (M f $ (M') f ® l) • (12) 

where |0) is one of the maximally entangled states 
{\GHZ) , |W)}. In this equation dHJ), the filters M, M' 
and the noise S act on different subsystems. This allows 
us to apply the evolution equation for bipartite concur- 
rence (|9|) to a "bi-partite" split 12 1 3 of the three-qubit 
system. We therefore obtain 

C 12 ^[(l®l®S)\x)(x\} = 

C 12 I 3 [(1®1®5)|0)(0|]C 12 I 3 [| X )], (13) 

while similar relations can be obtained for the "bi- 
partite" concurrences C 13 ' 2 and C 23 ' 1 of the three-qubit 
system. Although the Eq. (fT3"|) has similar structure to 
the evolution equation for bipartite systems ©, they 
differ by the maximally entangled state |0) in their 
right-hand sides: the maximally entangled state |0) = 

Yli=i K) ® K) /V^2 of the bipartite system is to be sub- 
stituted in Eq. (j9)), while one of the maximally entangled 
states (fTU|) - (fTTj) should be used in the right hand side of 

Eq. CED. 

Because of the symmetry of the maximally entangled 
states (fTU |) -(fTT |) with regard to the qubits permutation 
we have a relation 

C 12|3 [(1®1®S)|0)(0|] = C 13|2 [(1®S®1)|0)(0|] 
= C 23 I 1 [(5®1®1)|0)(0|], (14) 

where \<j>) = {\GHZ) , \W)}. From Eqns. USD and ft® 
it follows that for an arbitrary pure three-qubit state \x) 
the evolution of the bipartite concurrence is independent 
on a bipartite cut of the three-qubit system, i.e 

C 12 ^[(l®l®S)\x)(x\] = c 13 \ 2 {(i®s®i)\x)(x\] 
= C 23 ^[(S®l®l)\x)(x\}- (15) 

Substituting the evolution equation (|13[) into definition 
of the lower bound ^ and taking into account relation 
(|14p . we finally obtain an evolution equation of the lower 
bound for three-qubit concurrence 

r 3 [(l®l®<S) |x) (X|] - r 3 [(l®l®<S) |0) (0|]r 3 [|x)] , (16) 

where t 3 [..] is defined in Eq. ([5]). The entanglement dy- 
namics of an arbitrary pure state \x) of a three-qubit 
system, when one of its qubits undergoes the action of 
an arbitrary noisy channel 5, is subjected to the dy- 
namics of one of the maximally entangled states |0) = 



{\GHZ) , IT 7 ^)}. The choice between the maximally en- 
tangled states should be done after determining the en- 
tanglement class of the given state \x) following the pro- 
cedure in Ref. [24| and briefly discussed above. We note, 
that due to Eq. (|15[) the entanglement dynamics of a pure 
three-qubit state |x) is independent of which of the qubits 
is affected by the noise. In fact, this equation (|T5|) signifi- 
cantly simplifies the calculation of the lower bound ([5]). It 
is sufficient to compute just one bipartite concurrence in 
definition ([5]) of the lower bound, for example C 12 ' 3 [|x)], 
while the bipartite concurrences C 13 ' 2 [|x)] and C 23 ' 1 [|x)] 
are equal to it due to Eq. (|15p . 

D. Remarks on the evolution of pure iV-qubit and 
mixed three-qubit states 

It is desirable, of course, to generalize the evolution 
equation (|T6"| of the lower bound to the three-qubit con- 
currence also to ./V-qubit states, if just one of the qubits is 
affected by a noisy channel S. In contrast to the classifi- 
cation of the three-qubit states, however, it is not known 
until now how many and which entanglement classes 
eventually exist for qubit systems with N > 4, while 
some classification is available for N = 4 [2^, [26[ . It is 
therefore not directly possible to generalize Eq. (fTB)) for 
arbitrary pure states of N qubits. Nevertheless, some 
entanglement classes are known also for general pure N- 
qubit states, such as the GHZ- and W-class. If a given 
(pure) N-qubit state |x) belongs to the GHZ- or W-class, 
the evolution equation (|16|) of the lower bound can be 
extended to 

TttKi 8 "- 1 ®^)!*) (x|] = 

r Ar [(l^- 1 ®5)|0)(0|]r N [|x)], (17) 

where |0) denotes the corresponding maximal entangled 
7V-qubit state [13 

\GKZ) N = -L(\0f N + \lf N ) , (18) 
|M/) Ar = -^(|10,...,0) + |01,...,0) + |00,...,l)) . (19) 

We can further analyze the lower bound © to the 
three-qubit concurrence in order to understand the en- 
tanglement evolution in those cases where one starts al- 
ready with an initially mixed state pq. Exploiting the 
convexity of the lower bound (|16|) as a valid entanglement 
measure, we have r 3 [ (1 ® 1 ® S) p } = Ts[J2i Pi(l ® 1 ® 

sWhWU] < EiftT3[(i®i®s)|v>>,(V>U. Makin § 

use of this inequality in Eq. (|16[) . we obtain 

r 3 [(l ® 1 ® S)p ] < t 3 [(1 ® 1 ® S) |0) (0|] r 3 [p ] (20) 

for the evolution of the lower bound and for an initially 
mixed state. As before, we assume here that just one of 
the qubits is affected by the noisy channel S. Moreover, 
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the inequality (f20|) can be generalized for local two- and 
three-sided channels, i.e. to cases in which two or even 
all three qubits are affected by some local noise. For 
example, for a local two-sided channel Si (g> S 2 <8> 1 = 
(Si ® 1 <g> 1) (1 <g> S 2 <8> 1) we find 

r 3 [(Si ® S 2 ® l)po] < t 3 [(Si ® 1 ® 1) |0) (0|] 

xt 3 [(18 5 2 ® 1)|0) (0|]r 3 [po]- (21) 

ft is this particular form of Eq. (f2~Tj) that gives rise to a 
sufficient criterion for finite-time disentanglement of ar- 
bitrary initial states being subjected to local multi-sided 
channels 0. 



III. DYNAMICS OF THE INITIALLY MIXED 
STATE 

The three evolution equations (fT6|) . (|20|) and (|2Tj) pro- 
vide us with a powerful tool in order to describe the 
time-dependent entanglement dynamics of a three-qubit 
system under local noise. Although the assumption was 
made that only one qubit is affected by a noisy channel 
[apart from Eq. (|21[) ]. nothing more need to be known ex- 
plicitly about the time evolution of the system's state. In 
the present section, we make use of Eq. (|20[) in analyzing 
the entanglement dynamics of the initially mixed state 
([1]) for a variety of particular single-qubit noise models. 
Using the definition (fTO ]l -([TT ]l of the |GHZ) and \W) in 
Eq. ([T]), we can compute the various parts of inequality 
([2TJ]) analytically when one qubit undergoes the action of 
noisy channels. 

Indeed, there are several reasons for studying the en- 
tanglement evolution of the mixed state ([1]). For this 
state, first of all, an analytical expression is known for 
the convex roof to the concurrence (|3]) [l4[ . This enables 
one to compare the time-dependent lower bound from 
the evolution equation (|2U)) with the behavior of the con- 
vex roof as deduced from the state dynamics under the 
influence of a certain noise model. Second, the mixed 
state density matrix ((T|) has simply rank two. As we 
have recently shown in particular examples, for rank- 2 
density matrices the lower bound © coincides with the 
convex roof Third, it will be quite easy to compare 
also the speed of 'disentanglement' for an (initially) pure 
GHZ state p(t,p = 1) and the pure W state p(t,p = 0) 
under decoherence. We also note that the values for the 
lower bound for these pure states are related to each other 
through the ratio 

73(p(t = 0,p=l)) = M\GHZ}) = J_ 

Mp(t = o, P = oj) t 3 (\w)) 2V2' [ ' 

a result that was obtained in Ref. Q by means of a lower 
bound to the concurrence, different from definition ©. 

In the evolution equation (|2U|) we need to obtain the 
final state of the three-qubit system pg n under the influ- 
ence of a local noise S, i.e. p nn = (1 ® 1 ® S) pi n i, where 
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FIG. 1: (Color online) Evolution of the lower bound Tz(p) 
for initially mixed state {TJ if affected by the phase damping 
channel. While the blue surface shows the lhs of inequality 
(|20p . the red lines represents its rhs. 
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FIG. 2: (Color online) The same as in Fig. 1 but for the 
amplitude damping channel. 



Pini = {pa, \4>) (011- Following quantum operation for- 
malism [l[ , the final state can be obtained with the help 
of (Kraus) operators [28| as 

i 

and the condition J^. K[ Ki < I is fulfilled. The Kraus 
operators are well known for various single-qubit noise 
models; in the next section, we consider two frequently 
used noises, namely the phase and the amplitude damp- 
ing channels which are associated with relevant physical 
processes. 

A. Phase and amplitude damping channels 

Let us start the discussion with a three-qubit system 
which is prepared initially in the state ([1]) and for which 
just one qubits undergoes the action of the phase damp- 
ing channel. A phase damping describes for instance a 
diffusive scattering interaction of the qubit with its en- 
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vironment and is known to result into a loss of phase 
coherence information A possible representation of 
the phase damping in terms of time-dependent (Kraus) 
operators is given by [l| 

<=(%";), *r=( vT T !rr !!).(24) 

where T denotes a coupling constant. For this noise 
model, Fig. [T] displays the time-dependent evolution of 
the lower bound T3 for the initially prepared state (p} for 
different parameters p of the mixed state and at differ- 
ent times of the system-channel coupling. In this figure, 
the blue surface displays the left-hand side (lhs) of the 
inequality ([2"U]l . while the red lines shows the correspond- 
ing right-hand side. For all parameters < p < 1 of the 
mixed state, the lower bound T3 decays exponentially and 
vanishes only asymptotically for t — > 00. For the phase 
damping channel, moreover, the lhs and rhs of (|20|) are al- 
ways equal for an arbitrary parameter p and for all times 
t. 

If the same system is affected by a (local) ampli- 
tude damping channel, which describes the dissipativc 
coupling of a qubit to a thermal reservoir in the zero- 
temperature limit the operator elements in Eq. (|23|) 
are given by 

*!*-(; A). *?-{i*=r).m 

For such an amplitude damping, the time evolution for 
the lower bound T3 differs from the corresponding dynam- 
ics in the phase damping channel as seen from Fig. [2] In 
this amplitude damping model, in particular, the (blue) 
surface for the lhs of (f20f differs significantly from the 
rhs of this inequality for some values of the parameter p. 

B. Sudden death of three-qubit entanglement 

Based on violation of additivity in the case of entan- 
glement decay in a multi-qubit system coupled to two 
independent weak noises, entanglement sudden death Q 
reveals a practically important aspect of time-dependent 
entanglement evolution. It is important to verify whether 
such a phenomenon can be predicted with the suggested 
evolution equation of the lower bound (|20|) . Let a three- 
qubit system be prepared in the state (Q]) and one of the 
qubits is subjected to the generalized amplitude damping 
channel. This noisy channel can be viewed as a 'superpo- 
sition' of two independent amplitude damping channels 
acting on a qubit and can be expressed by four Kraus 
operators Kf ad = \Kf, Kf & = \Kf and 

K T & = \ ( " 0" J ) , Kf d = \ (Kff , (26) 

where Kf d and Kf d are defined by Eq. ([231) . The evolu- 
tion of the lower bound T3 for the three-qubit state ([1]) 




1.0 



FIG. 3: (Color online) The same as in Fig. 1 but for the 
generalized amplitude damping channel. The sudden death 
of entanglement is clearly seen for all < p < 1. 



is shown in Fig. [3l Although the lhs and the rhs of the 
inequality (|20[) differ for some parameters p of the initial 
state p{p) , they both vanish in a finite time for all values 
P- 



IV. CONCLUSION 

Unlike the evolution equations of concurrence for two 
qubits Q and for a bipartite system [{|, we have pre- 
sented an evolution equation of the lower bound for 
multi-qubit concurrence (|16p for an initially pure state 
of the system and when just one qubit is affected by lo- 
cal noise. This evolution equation (fTl)]) has been also ex- 
tended to the cases of initially mixed states (f2"U)) and the 
action of many-sided noisy channels (|2ip . In addition, 
the evolution equation of the lower bound for initially 
mixed states (|2TJ|) has been employed especially to show 
entanglement dynamics of the mixed state (HJ) when just 
one qubit is affected by the phase, the amplitude or the 
generalized amplitude damping channel [cf. Figures [I] [3]. 

Of course, the lower bound ((4]) is only an approxima- 
tion to the convex roof for multi-qubit concurrence §3§ 
for which an analytically computable expression is cur- 
rently unavailable. Nevertheless, we have recently shown 
that this lower bound coincided with the convex roof for 
some scenarios of entanglement dynamics p|. A more 
detailed analysis of the accuracy of the lower bound ap- 
proximation is presently under work. 
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